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ABSTRACT 

We investigate the possible effects of the supernova ejecta hitting the companion star in iPTF 13bvn, 
focusing on the observable features when it becomes visible. iPTF 13bvn is a type lb supernova that 
may become the first case that its progenitor is identified as a binary by near future observations. 
According to calculations by Bersten et al. (2014), the progenitor should have a mass ~ 3.5 Mq to 
reproduce the supernova light curve, and such compact stars could only be produced via binary 
evolution. This is one of the reasons that we expect the progenitor to be a binary, but it should be 
confirmed by observing the remaining companion after the supernova. Their evolutionary calculations 
suggest that the companion star will be an overluminous OB star at the moment of supernova. With 
a combination of hydrodynamical and evolutionary simulations, we find that the secondary star will 
be heated by the supernova ejecta and expand to have larger luminosities and lower surface effective 
temperatures. The star will look rather like a red super giant, and this should be taken into account 
when searching for the companion star in the supernova ejecta in future observations. 

Subject headings: binaries: close — stars: evolution - supernovae: individual (iPTF 13bvn) 


1. INTRODUCTION 

Type lb supernovae (SNe) are one of the hydrogen de¬ 
ficient subtypes of core-collapse supernovae, which are 
the final fates of massive ( M > 8Mq) stars (Filippenko 
1997). Their lack of hydrogen indicates that they ex¬ 
plode from stripped-envelope progenitors. Two major 
scenarios have been proposed to explain the removal 
of the hydrogen layers. One is by strong stellar winds 
for stars with zero age main sequence (ZAMS) masses 
Mzams > 25 Mq (Maeder 1981). Such stars are called 
Wolf-Rayet stars, and are known that they sometimes 
shed their entire hydrogen envelope with the wind. The 
other possible scenario is mass transfer in close binaries 
(Podsiadlowski et al. 1992). Outer layers of the more 
evolved star in a binary can be transferred to its compan¬ 
ion via Roche lobe overflow, or removed dynamically in 
common envelope phases. Debates continue about which 
is the more likely scenario. 

iPTF 13bvn was a recent SN of this particular type, 
first identified by the intermediate Palomar Transient 
Factory (Law et al. 2009) on June 16.238 UT 2013. 
Its host galaxy is NGC 5806, which is at a distance 
~ 21Mpc. Pre-explosion images of Hubble Space Tele¬ 
scope (HST) showed a candidate progenitor within a 2a 
error of the SN site (Cao et al. 2013). It is not settled 
whether we were looking at the progenitor star itself or a 
combined flux of a binary, but nevertheless we can place 
strong constraints on the exploding star. 

Groh et al. (2013) first proposed that a Wolf-Rayet 
star with a ZAMS mass of 31 — 35 Mq could be a possible 
progenitor for iPTF 13bvn. This result was based on the 
absolute magnitude of the source in the pre-explosion 
image, and the strict upper limit on the radius of the 
progenitor (< 5 Rq) due to early detection. In their 
case, the pre-SN progenitor mass was ~ 10 Mq. This 
possibility was ruled out by observations of the later 
phases of the bolometric light curve (Bersten et al. 2014; 
Fremling et al. 2014). Detailed hydrodynamical simula¬ 


tions of SNe with different progenitors showed that the 
observed light curve cannot be reproduced by a star with 
M > 8 Mq. Instead, the preferred model was a He star 
with M ss 3.5 Mq which is difficult to produce by assum¬ 
ing only stellar winds. The limit on the stellar radius was 
also extended up to < 150-Rq by an additional set of sim¬ 
ulations with extended thin envelopes by Bersten et al. 
(2014) (hereafter B+14). Their results strongly support 
the binary evolution scenario, and they also showed a 
possible evolutionary path that can produce such a low 
mass progenitor star and match the pre-SN HST obser¬ 
vations consistently. The SN is estimated to fade below 
the brightness of the pre-SN primary star at about three 
years from explosion. A fainter secondary star may be 
observed shortly afterwards. If a companion star is really 
discovered, it will be the first case that the progenitor for 
a type lb SN is confirmed to be a binary. 

According to the binary evolution calculations by 
B+14, the secondary star is predicted to be an overlu- 
minous OB-type star. The combined optical flux of the 
primary and the secondary was consistent with the opti¬ 
cal flux from the pre-SN image. However, the SN ejecta 
may change the observable features of the companion af¬ 
ter the explosion by stripping mass, or heating the star 
to make it swell up (Hirai et al. 2014). This effect was 
not included in previous predictions although it may sig¬ 
nificantly alter its appearance. 

In this paper we investigate the possible effects of SN 
ejecta hitting the companion in iPTF 13bvn, and how it 
may affect the observational features. Such events have 
been thoroughly studied for the single-degenerate sce¬ 
nario of type la SNe (Marietta et al. 2000; Podsiadlowski 
2003; Pan et al. 2013; Liu et al. 2013; Shappee et al. 
2013), but not so much for massive binaries, and rela¬ 
tively wide binaries, which do not fit in the picture of the 
single-degenerate scenario. Here we first model binary 
systems that fit within the constraints placed by B+14 
via stellar evolution calculations, and then carry out hy- 
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drodynamical simulations of the SN ejecta colliding with 
the companion star. To obtain the stellar structures at 
the time it becomes visible in the ejecta, we then perform 
additional stellar evolution simulations with extra heat 
distributed in the outskirts of the star. Parameters for 
the artificial heating were evaluated from the hydrody- 
namical simulations. 

This paper is structured as follows: Our choice of 
model parameters and numerical method is explained in 
Section 2. The results of our evolution and hydrodynami- 
cal simulations are explained in Section 3, and we discuss 
our results in Section 4. We summarize our conclusions 
in Section 5. 

2 . MODELS AND NUMERICAL METHOD 

B+14 revealed that the progenitor mass of iPTF 
13bvn should be in the range of 3-5M 0 , presumably 
~ 3.5M 0 . By detailed binary evolution calculations, 
they also placed constraints on the initial configuration 
of the binary system that can produce such stripped 
envelope progenitors in this mass range. According 
to their estimates, the initial primary star (progenitor) 
mass should be 15 < M[ < 25 Mq, and the mass ra¬ 
tio 0.8 < M\/M\ < 0.95. Note that the lower bound 
of the mass ratio is not physically motivated, but set 
to avoid common envelope phases since they could not 
handle common envelope evolution with their stellar evo¬ 
lution code. These constraints overlap with the results 
by Eldridge et al. (2015), where the permitted primary 
mass was IOMq < M{ < 20Mq. The overlapping range 
may be most plausible. They further estimated the pos¬ 
sible range of the final secondary mass, being 23 < M/ < 
45 Mq if conservative mass transfer is implied. This range 
is expanded to 18 < Ml < 45 Mq if mass accretion effi¬ 
ciencies are taken in the range 0.5 < /? < 1 where /3 is 
the fraction of transfered matter that is accreted onto the 
secondary star (Benvenuto et al. 2013). The luminos¬ 
ity of the secondary star ended up within the constraint 
placed by the pre-explosion image of iPTF 13bvn. 

In this paper, we first attempt to construct a pre-SN bi¬ 
nary model that fits within the observationally permitted 
range. We utilize the open source stellar evolution code 
MESA (version 7184; Paxton et al. 2011, 2013) to model 
each star in the binary. Binary evolution is treated by 
evolving two spherically symmetric one dimensional stars 
with a mass transfer rate applied by the “Ritter formal¬ 
ism” (Ritter 1988). Orbital elements are also evolved 
consistently. We assume no rotation, no magnetic fields, 
and circular orbits for all models to maintain simplicity. 

It is known that even for identical initial conditions, 
different stellar evolution codes give different results 
(Jones et al. 2014; Sukhbold & Woosley 2014). Since the 
code employed in B+14 (a code by Benvenuto et al. 2013; 
Benvenuto & de Vito 2003) was different from ours, we 
could not produce the same binary model as the example 
case in B+14 by taking the same initial conditions. Here 
we aim to reproduce similar pre-SN configurations as the 
one shown in B+14, so we take slightly different initial 
parameters. We also construct several other binary mod¬ 
els for comparison. The selected binary parameters are 
listed in table 1 , along with the final binary parameters at 
the moment of SN. For the non-conservative mass trans¬ 
fer models (with /? = 0.5), the angular momentum loss 


Table 1 

Binary parameters of the stellar evolution calculations with MESA. 


Model 3, 

Age 

(Myr) 

Mi 

(M 0 ) 

m 2 

(M 0 ) 

Ri 

(R&) 

R 2 

(Rg) 

P 

(days) 

a 

(Rg) 

al.O 

0 

19.0 

18.0 

5.65 

5.47 

2.45 

25.5 


10.4 

3.48 

33.5 

44.7 

7.72 

61.7 

219 

a0.5 

0 

19.0 

18.0 

5.65 

5.47 

2.45 

25.5 


10.2 

3.72 

25.6 

47.8 

10.4 

62.0 

203 

bl.O 

0 

18.0 

17.0 

5.47 

5.30 

3.50 

31.7 


9.98 

4.15 

30.9 

15.4 

8.02 

47.9 

181 

b0.5 

0 

18.0 

17.0 

5.47 

5.30 

3.50 

31.7 


9.94 

4.23 

23.9 

10.6 

7.76 

54.4 

184 

cl.O 

0 

15.0 

14.0 

4.93 

4.74 

3.30 

28.7 


12.5 

3.02 

26.0 

31.8 

7.01 

63.4 

206 

c0.5 

0 

15.0 

14.0 

4.93 

4.74 

3.30 

28.7 


12.5 

3.06 

20.0 

30.4 

6.57 

73.3 

210 


Note. — All models were assumed to have circular orbits. 


a Numbers in the model names indicate values of the accretion 
efficiency parameter [3 (Benvenuto et al. 2013) applied in our evo¬ 
lution calculations. 

parameter is given as a = 1 (Benvenuto & de Vito 2003). 
All models were calculated until itthe primary starts pro¬ 
ducing iron at the centre. The primary stars are assumed 
to explode within hours after iron production, and binary 
parameters are unchanged until explosion. 

To the pre-SN binary model, we then apply the same 
method as in Hirai et al. (2014) to simulate the SN ejecta 
hitting the companion. A two step strategy is employed; 
1 ) artificially explode the primary star on a one dimen¬ 
sional spherical grid, 2) simulate the effect of SN ejecta 
hitting the companion star on a two dimensional axisym- 
metrical grid. The hydrodynamical code “yamazakura” 
(Sawai et al. 2013) is used for all hydrodynamical simu¬ 
lations throughout this paper. It is a mesh-based central 
scheme code with an iterative Poisson solver to solve self¬ 
gravity. An ideal-gas equation of state was applied with 
an adiabatic index 7 = 5/3. 

For the first step, E^md + E exp of internal energy is 
added to the inner few meshes of the primary star model 
placed on a one dimensional computational domain to 
initiate an explosion like a thermal bomb (Young & Fryer 
2007), where E^nd is the gravitational binding energy of 
the star and E exp is the explosion energy. In order to 
leave a residual neutron star after the explosion, we cut 
the central 1.4 Mq of the star out of the computational 
region, and set a reflective inner boundary condition. In 
this way, all the energy applied will be placed just above 
the neutron star and it will initiate a shock wave that 
propagates outwards through and out of the star. The 
explosion energy was taken as E exp = 8 x 10 5O erg from 
the estimate by B+14, which was obtained to fit the peak 
of the bolometric light curve and photospheric velocity 
evolution (See also Srivastav et al. 2014, for discussions 
on the explosion energy). A dilute circumstellar matter 
is placed around the star due to numerical reasons, with 
a low density so as not to disturb the propagation of 
the SN ejecta. An outgoing outer boundary condition is 
employed so the ejecta can flow out freely. Ejecta profiles 
are sampled at a point far from the stellar surface (50 
times the stellar radius), and it is checked that it follows 
a homologous expansion by comparing the time evolution 
at different points. 
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In the second step we place the secondary star model 
onto the origin of an axisymmetric 600(r) x 180(0) spher¬ 
ical grid. The axis is taken along the line connecting the 
centre of both stars. Axisymmetry is justified due to the 
short timescale of SN ejecta flowing past the star (~ 1 
day) compared to the orbital period (~ 60 days). Un¬ 
like Hirai et al. (2014), we do not leave out the central 
portion of the star since the companion star model is 
not so centrally concentrated. The density scale height 
is resolved with at least 10 meshes, with a total ~ 150 
radial gridpoints inside the star. Mesh sizes are then in¬ 
creased monotonically outside the star, until the radius 
of the computational region reaches up to ~ 90% of the 
binary separation. Yet again we place a low density at¬ 
mosphere around the star that is dilute enough so that 
it has negligible mass. Data from the first step is used to 
place the SN ejecta close to the companion as an initial 
condition. For the outer boundary condition, the ejecta 
data is extrapolated as a Dirichlet boundary on the side 
facing the primary, whereas a free boundary is applied 
on the opposite side. 

Each mesh is marked as bound or unbound using the 
“Bernoulli criterion”, in which matter is bound when 

\v 2 +£+-+(/)< 0 (1) 

2 p 

where e is the specific internal energy and p is pressure, 
p is density and 4> is the gravitational potential. We 
integrate over the bound region to evaluate the mass and 
position of the centre of mass of the remaining star. In 
order to see whether matter from the ejecta can mix with 
the original stellar matter, we also placed tracer particles 
that just follows the fluid motion. Each particle carries 
information of its mass and origin (stellar or ejecta), and 
is evaluated whether it is bound or not at each step. 

3. RESULTS 

3.1. Stellar evolution 

In this section we show results of our stellar evolution 
calculations. Figure 1 shows the evolutionary track on 
the HR diagram of both stellar components for models 
al.O and a0.5. They roughly resemble the track of the 
model presented in B+14. The overall evolution of the 
primary does not change much for different values of /3, 
nor does the period of the system. On the other hand, 
the secondary mass and luminosity depends strongly on 
the value of /3. However, the optical flux is dominated by 
the primary 1 , so it is difficult to constrain the secondary 
luminosity from the HST pre-SN image. All other models 
followed similar evolutionary tracks on the HR diagram. 

Every primary star model was deficient of hydrogen 
at the point of SN, which is a necessary condition to 
produce a type lb SN. Final luminosities of the primary 
stars range 4.5 < L,{/Lq < 4.9, which roughly fits in 
the allowed range to match the HST pre-SN flux. This 
is also the case for the companions, with luminosities 
4-9 < L^/Lq < 5.6. From the above facts, we assume 
that all of our models can be progenitor candidates of 
iPTF 13bvn. Since the progenitor mass is most likely 
ss 3.5Mg, we will take model al.O to be our reference 

1 Even though the secondary star has a higher luminosity, the 
primary has a stronger optical flux due to the lower temperature. 



log (T eff !K) 


Figure 1. Evolutionary tracks of the stellar components in models 
al.O (black lines) and a0.5 (red lines). Solid lines show results of 
the primary stars, and dashed lines the secondary stars. The star 
signs label the position of the primary star at the point of SN. 



Figure 2. Density profiles of the companion star at the point of 
SN of the primary. Model names correspond to the names in table 
1. 

model. The b series have primary masses close to the 
lower limit of the observational constraint and the c series 
close to the upper limit. 

Density-radius profiles for each companion star model 
are compared in Figure 2. All models have similar struc¬ 
tures, with steep density gradients near the surface and 
slightly different radii according to their accretion histo¬ 
ries. 

Obviously, these are not the only models which can re¬ 
produce the observational characteristics of iPTF 13bvn. 
Common envelope phases may have occured to produce 
the compact progenitor, or the orbit may have been ec¬ 
centric to induce periodic mass loss. Single star mod¬ 
els cannot be excluded if we consider stellar rotation 
(Eldridge et al. 2015). Our stellar evolution calculations 
do not include these effects, and it is out of the scope of 
this paper to consider every possible scenario. 

3.2. Supernova 

Figure 3 shows snapshots of our hydrodynamical sim¬ 
ulations of the collision of SN ejecta and the companion 
star in our reference model. Note that the scales are 
different in the left and right panels. The whole compu¬ 
tational region is displayed in the right panels whereas we 
show a close-up view in the left panels. Panel a shows 
our initial condition, where the ejecta is just about to 
reach the companion surface. As the ejecta hits the sur¬ 
face of the companion, a forward shock is driven into the 
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star, and a reverse bow shock is formed in front of the 
star which can be seen as density discontinuities in Panel 
b. The forward shock heats the outskirts of the star as 
it proceeds, but the propagation soon stops at a certain 
point because of the pressure gradient in the star. After 
the bulk of the ejecta flows past, the bow shock and the 
heated matter starts to expand outwards as can be seen 
in the lower half of Panel c. We can check whether the 
stellar matter is expanding, or it is just the bound region 
that is expanding by looking at the tracer particle distri¬ 
butions. Figure 4 shows the distributions of the tracer 
particles at times corresponding to Panels b and c in Fig¬ 
ure 3. The upper half is colour coded according to the 
origin of each particle, red for ejecta matter that comes 
from the primary, and blue for secondary star matter. 
The lower half is colour coded whether it is bound (light 
blue) or not (grey). We can see that some of the parti¬ 
cles from the stellar surface is stripped away along the 
axis, but the mass is very small. If there are any parti¬ 
cles that are red in the upper half and light blue in the 
lower half, it means that ejecta matter has accreted onto 
the companion. However, most of the light blue parti¬ 
cles in the lower half are blue in the upper half, which 
indicates that the enlargement of the bound region is not 
due to ejecta matter being accreted to the stellar surface, 
but the shocked stellar matter expanding because of the 
heat. A small amount of ejecta matter is mixed into the 
stellar matter, but their masses are extremely small. The 
bound region keeps on expanding in a spherical manner, 
eventually reaching the outer boundary of our simula¬ 
tion. Panel d in Figure 3 shows a snapshot at these later 
stages, where the bound region is almost spherical. 

In Figure 5 we show the amount of unbound mat¬ 
ter (M u b), and the displacement of the centre of mass 
( A xcom) as a function of time since SN. As we can see 
from the upper panel, hardly any matter was stripped 
off nor accreted onto the star. The absolute values are 
comparable to the errors arising from the limitation in 
numerical resolution. This is consistent with previous 
works (Hirai et al. 2014; Marietta et al. 2000; Pan et al. 
2010 ), where the unbound mass depends on the orbital 
separation as oc a m “ 6 , where a is the separation. The 
exponent m u b depends on the stellar structure, but sim¬ 
ulations suggest values —4.5 < m u b < —3. Since our 
models all have relatively wide separations, it is natu¬ 
ral that the amount of stripped mass was so small. The 
lower panel shows the motion of the centre of mass, which 
reaches a constant velocity in the direction opposite to 
the exploding star in about a few days. According to sim¬ 
ple estimates (Cheng 1974; Wheeler et al. 1975), this so- 
called “rocket effect” or “kick” velocity 2 should be very 
small in our present model, due to the wide separation 
with respect to the stellar radius. In our simulation, the 
kick velocity reached up to 1.75 x 10 4 km h _1 , which was 
~ 30% of the orbital velodity. This relatively large kick 
velocity may help destroy the binary or at least alter the 
eccentricity of the post-SN system. For model a0.5, the 
kick velocity was ~ 5% of the orbital velocity, which is 
almost negligible. All other models showed qualitatively 
similar results. 

In Table 2 we list the amount of unbound mass and 

2 The kick velocity mentioned here is different from the natal 
kick imparted to the central neutron star in core collapse SNe. 


the kick velocity along with the final temperature and 
luminosity which is explained in Section 4.1. It is clear 
that the unbound masses were negligible in every model. 
The kick velocity was also negligible compared to the or¬ 
bital velocities. There seems to be no strong correlation 
between the binary parameters to the unbound mass or 
kick velocity. This is because the stellar structure and 
SN ejecta profiles are different for each model, and the 
absolute values are all comparable to the numerical res¬ 
olution. 


Table 2 

Final states of the secondary star. 


Model 3, 

Mub 

(M 0 ) 

v kick 

(km h 1 ) 

^orb 

(km h -1 ) 

log T eff c 
(K) 

logL° 

(L G ) 

al.O 

2.1 x 10 -4 

1.8 X 10 4 

6.1 x 10 4 

3.82 

5.51 

a0.5 

4.7 x 10~ 4 

4.3 x 10 3 

7.6 x 10 4 

3.64 

5.65 

bl.O 

1.1 x 10“ 4 

5.6 x 10 3 

8.2 x 10 4 

4.01 

5.85 

b0.5 

3.1 x 10“ 4 

9.1 x 10 3 

9.3 x 10 4 

3.66 

5.58 

cl.O 

0.9 x 10 -4 

3.3 x 10 3 

6.1 x 10 4 

3.61 

5.57 

c0.5 

2.2 x 10“ 4 

5.7 x 10 3 

6.9 x 10 4 

4.22 

4.98 


“Names of models correspond to those in Table 1. 
b The orbital velocity prior to SN. 

c T e ff and L are values estimated three years after SN. 


4. DISCUSSIONS 

The dynamical effects of the ejecta hitting the compan¬ 
ion star, such as mass stripping and momentum trans¬ 
fer seem to be limited. But the expansion of the outer 
layers of the star may become important when we ob¬ 
serve the companion after the SN ejecta becomes faint 
enough. Such expansions have also been found in previ¬ 
ous simulations (Hirai et al. 2014; Liu et al. 2013). Here 
we explore the observable features of the companion star 
when it becomes visible. 

4.1. Reddening due to SN Heating 

In Figure 6 we illustrate the evolution of the total en¬ 
ergy of the star, which we calculate as 

Etot = [ {(f> + e + v 2 )pdV (2) 

Jv b 

where is the bound region. The SN ejecta imparts 
energy to the star via shock heating, and the star adjusts 
its structure to restore mechanical equilibrium. After a 
few days, the total energy of the star reached a constant 
value. We consider that the heating process has finished 
here, and assume that the difference between the initial 
and final energy is the amount of energy injected to the 
star by the SN ejecta. 

In our hydrodynamical simulations, the expansion of 
the star did not cease even at the end of our simula¬ 
tions (a few days) where the bound region reaches the 
outer boundary. The numerical cost is extremely expen¬ 
sive to wait for the star to recover equilibrium. Hence 
we use the stellar evolution code instead of the hydro- 
dynamical code to estimate the state of the star for 100 
years after SN. We believe that this is justified because 
the expansion was almost spherical in the hydrodynam¬ 
ical simulations. This is a similar procedure to previous 
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Figure 3. Density plots from the hydrodynamical simulations of the collision of SN ejecta and the companion star. Each snapshot is 
labelled with the time ellapsed since SN. Only the bound matter is coloured in the lower half of each panel. The SN ejecta is flowing in 
from the left side of each panel. The left panels have a radius of 1.3 X 10 'km, and right panels 1.3 X 10 8 km. 


studies that predicted the post-SN appearance of the 
companion in type la SN progenitors in the single de¬ 
generate scenario (Podsiadlowski 2003; Pan et al. 2012; 
Shappee et al. 2013). 

Starting from the companion model at the pre-SN 
stage of the primary , we artificially input energy into 
the outer layers as 


Cheat (m) 


A E heat ( (m — p) 2 \ 

t heat V2^ e ^ \ 2a^ ) 


( 3 ) 


where Cheat is the artificial specific internal energy in¬ 
jection rate as a function of mass coordinate, A Eh ea t is 
the total energy input, theat is the duration of heating, 
and a = (M/ - m{r in ))/ 6, pt = (Af/ + m(r in ))/2. r in 
is the radial coordinate at the inner edge of the heating 
layer, where the shock could not proceed further inwards 
in the hydrodynamical simulations. Figure 7 illustrates 
the time evolution of the distribution of specific internal 
energy along the axis. The black line shows the initial 
condition, where we can see the ejecta approaching the 
stellar surface (~ 5.3 x 10 6 km) from the left. The forward 
and reverse shocks can be seen as a steep wall in the red 
line, and weakens its strength as it proceeds into the star 
which can be seen in the blue line. The forward shock 


propagates up to ~ 3 x 10 6 km from the centre at most. 
We can see from the green line, the internal energy dis¬ 
tribution near the end of our simulation, that the stellar 
matter interior to this point is almost unaffected by the 
shock. We therefore take rt n = 3 x 10 6 km as the inner 
edge of the heating layer. AEh ea t and theat are also set 
to the values taken from the hydrodynamical simulations 
( AEh ea t = 1 x 10 5O erg, t hea t = 100hrs). Our strategy 
of extra energy deposition is similar to Shappee et al. 
(2013) and Podsiadlowski (2003). The main difference 
is that we ignore the effects of mass stripping, since the 
amount of stripped mass in our hydrodynamical simula¬ 
tions were negligible (cf. tabic 2). Due to this assump¬ 
tion, all the energy from the SN ejecta is used to puff up 
the star, not to ablate away some surface matter. 

The shock heating of the outer layers by the SN ejecta 
modifies the stellar structure, pushing out the outer lay¬ 
ers like it did in the hydrodynamical simulations. In Fig¬ 
ure 8 we show the evolutionary track of the artificially 
heated companion star on the HR diagram for 100 years 
after SN. Focusing on the fiducial model (black line), it 
can be seen that the temperature decreases almost mono- 
tonically. During the heating phase, temperature and lu¬ 
minosity both drops but the luminosity eventually starts 
to increase due to the growing expansion speed. As soon 
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Figure 4. Distributions of tracer particles at various times. In the upper half of each panel, each particle is coloured red or blue if they 
originate from the ejecta or star respectively. In the lower half, particles are coloured light blue or grey depending on whethere they are 
bound or not. Particles that are red in the upper half and light blue in the lower half are ejecta matter that has mixed into the stellar 
matter. The radii of the circles are 1.15 X 10 8 km. 



Figure 5. Time evolution of the fraction of unbound mass (upper 
panel) and the displacement of the centre of mass (lower panel). 
Red lines show results for model al.O, and blue lines for a0.5. Black 
slopes in the lower panel show the approximate velocities of the 
centre of mass. 

as we switch off the artificial heating, the temperature 
and luminosity shows an almost discrete drop because of 
the disappearance of heat 3 . In the following few dynam¬ 
ical times, the radius and luminosity increases to restore 
hydrostatic equilibrium. We do not expect these large 
fluctuations in the early phase to be realistic, since it 

3 The artificial heat is applied all the way up to the surface, 
contributing to the surface temperature. 



Figure 6. Evolution of the total energy of the star. Total energy 
is evaluated by Eq.2 at every time step. The difference between 
the initial and final energies are considered to be the heat imparted 
by the SN ejecta A E heat . 

is highly affected by our methodology. Thus we only 
pay attention on the evolution after several dynamical 
times. After a year, the effective temperature will drop 
by an order of magnitude, and the luminosity slightly 
increases. It was roughly estimated in B+14 that three 
years is the time it takes for the SN light curve to decline 
below the brightness of the original progenitor. The time 
it takes for the companion to become visible depends on 
the declination of the SN light curve and the band used 
for the observation, but we expect that it is not so far 
from three years, maybe slightly earlier. At that time, it 
will be somewhere around the middle star sign in Figure 
8 , somewhat like a red supergiant. Even after it becomes 
visible, we expect the star to continue increasing its lu¬ 
minosity for a few decades. 

Although we have carefully chosen our parameters ?y n , 
AEheat, and theat to match our hydrodynamical simu¬ 
lations, there are still some uncertainties left. Results 
for the same simulation with different parameters are 
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Figure 7. Time evolution of the distribution of specific internal 
energy along the axis. The direction of the axis is taken away from 
the exploding star. 

also shown in Figure 8. The final temperature strongly 
depends on ri „, where larger values of rj n led to lower 
temperatures and lower luminosity. This is because if ri n 
is larger, the extra energy is injected into a much smaller 
volume, so it needs to expand more to retain equilibirum. 
Final states are also sensitive to the value of A Eheat, 
which can be seen in Figure 8 where results for differ¬ 
ent A Eheat are plotted with different colours. Larger 
values of AEheat give lower temperatures and higher lu¬ 
minosities. We have also conducted the same simulation 
with different theat (black dotted line in Figure 8), but 
the differences were limited. As explained earlier, the 
early evolution is largely affected by our methodology, 
particularly on theat, but the differences become indistin¬ 
guishable within a few dynamical times after the heating 
phase. We are therefore sure that the long term evolu¬ 
tion will not depend on our choice of theat, and any ear¬ 
lier evolution (which is not observable anyway) should be 
studied by hydrodynamical simulations since it is highly 
unspherical. Regardless of the uncertainties in our as¬ 
sumptions, the secondary star will appear as a red star 
when it becomes visible with temperatures lower by a 
factor of ~2-10, off the main sequence. The parameter 
dependence of the luminosity is weak, just slightly in¬ 
creasing in some cases by less than a factor of ~2. 

Estimates of the final temperature and luminosity for 
all models are listed in Table 2. All other models have 
considerably lower temperatures compared to their orig¬ 
inal states. The expansion depends strongly on the stel¬ 
lar structure so it may be possible to distinguish among 
models from observations. 

4.2. Other Consequences 

The expansion of the star exceeded the original binary 
separation in some of our models, which means that it 
may engulf the newly born neutron star. For such cases, 
the pulsar wind emitted by the neutron star may heat or 
strip the loosely bound envelope and distort the shape. If 
the neutron star becomes embedded in the expanded en¬ 
velope, it might initiate a common envelope phase, pos¬ 
sibly leading to much tighter orbits or the formation of 
Thorne-Zytkow objects (Thorne & Zytkow 1977). It will 
be difficult to confirm whether the system has entered a 


Figure 8. Evolution of the secondary star after SN of the primary 
for model al.O. The original position is the left end of the lines, 
and the right ends are 100 years after SN. Star symbols mark the 
location at 1, 3 and 10 years since SN. Different colours show results 
with different AEheat ■> and the dashed and dot-dashed lines show 
results with different ri n . The black dotted line shows the evolution 
of the star with the same heating parameters as the black solid line, 
but with theat = 0.2 years. 

common envelope phase by observations, but if the red 
super giant-like companion is found and no neutron star 
is detected for the following couple of years, it may be 
the first ever evidence of an ongoing common envelope 
phase. 

5. CONCLUSION 

iPTF 13bvn is a candidate of the first case that the 
progenitor for a type lb SN is confirmed to be a bi¬ 
nary. There are several possible evolutionary paths that 
can both reproduce the SN light curve and the pre-SN 
photometry of iPTF 13bvn, all containing a compact 
(< 5M 0 ) He star primary and an overluminous OB sec¬ 
ondary star. We have studied the effects of SN ejecta 
colliding with this secondary OB star, and we found that 
the surface of the star may swell up with the heat, be¬ 
coming more red and luminous when it becomes visible 
in the SN remnant. The change in temperature and lumi¬ 
nosity strongly depends on the parameters A Eheat and 
ri n that we used in the post-SN evolution calculations. 
But in any case the secondary will be found off the main 
sequence, with lower effective temperatures and a slight 
increase in luminosity which also depends on the time we 
detect it. Effects of mass stripping and kicks are limited, 
and will not affect the further evolution of the binary. 

The radius of the secondary star becomes so large that 
it may engulf the primary-produced neutron star. In that 
case we might be able to observe the first ever ongoing 
common envelope phase as an absence of the neutron 
star in the vicinity. 

We note that other evolutionary paths are possible if 
we include stellar rotation, common envelope phases or 
eccentric orbits. Theories that can deal with all of these 
cases have not been established yet and it is much out 
of the scope of this paper (See Ivanova et al. 2013, for 
a review on common envelopes). If we concern these 
possibilities, there is a much wider range of possible pro¬ 
genitors so further studies are certainly warranted. 
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